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Abstract 

The calculation of the charge present in central region of the double barrier struc- 
ture at non-equilibrium conditions is discussed. We propose here a simple method 
to calculate non equilibrium Green's functions which allows consistent calculations 
of retarded and distribution functions. To illustrate the approach we calculate the 
charge on the quantum dot coupled via tunnel barriers to two external leads having 
different chemical potentials fj,L and /xr. The obtained results have been compared 
with other approaches existing in the literature. They all agree in the equilibrium 
situation and the departures grow with increasing the difference — Hr- 
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Recent advancements of nanotechnology allow the design and study of devices in 
which quantum phenomena play primary role IJ. Typical device consists of a small 
central region coupled via tunnel barriers to two external electrodes (leads). In such 
quantum-dot devices Coulomb blockade, resonant tunneling [[| and Kondo effects ||, 
have been observed. Theoretical description of electron transport phenomena in quantum 
dot often requires (self-consistent) calculations of the charge accumulated in the central 
region |j[ If the system is in an equilibrium state (i.e. the chemical potentials of the 
left and right leads are the same = Hr = the calculation of charge on the dot 
is easy. One finds the (local) density of states for the dot N^e) and integrates it with 
appropriate distribution function, which in equilibrium coincides with Fermi-Dirac one; 
f(e) = (ePie-ri + i)-i ; W ith p = l/k B T and gets 



n = J deN{e)f{e) . (1) 

— oo 

The easiest way to derive ([!]) is to find the (equilibrium) Green's function |7j of the dot 
and apply spectral theorem to get n. In non equilibrium situation (/iL 7^ Hr) one can not 
use the equilibrium technique. Physically the problem is much more complicated as the 
tunneling current flows across the system and there is no a priori simple prescription as 
how to calculate the charge on the central region of the device (dot). Some additional 
complications arise if time dependent voltage or other time dependent fields are applied 
to the system. Here, however, we shall consider non-equilibrium but steady state. At 
nonequilibrium one can think about at least two ways of generalizing Eq. ([!]). The first 
is to replace the chemical potential fi by suitably chosen average chemical potential ft. 
The other is to replace f(e) by the average distribution function f(e). The simplest 
two possibilities, which in fact have been used in literature || to discuss various issues 
connected with transport across the system at hand are the following: 

/i = ^— (2) 

or 

/(,)- /'(«) + /"M , (3) 

where /l/r — (e^ - ^/-^ + 1) _1 , is the equilibrium distribution function for the left/right 
lead. The immediate objection against the above proposals is that both give slightly 



different results, except at equilibrium and none of them is well justified. Equal weights 
assigned to Hl(r) or fL{R){^) may possibly be also justified for symmetrically coupled 
systems only 

Below we shall present an approach which allows the correct calculation of the charge 
n in non equilibrium systems. To make the presentation of it as simple as possible we 
shall resort to the specific model describing the quantum dot coupled to two external 
leads. We use here the Anderson-Hubbard model 

H = Y,( £ ^ ~ ^\)C\kaC\ka + E d Y,dtd a + Un 1 n l + £(^4^ + h.c.) . (4) 

Afccr a \ka 

Here A = R,L denote the right (R) or left (L) lead in the system. The parameters have 
the following meaning: c\ ka (c\k a ) denote creation (annihilation) operator for a conduction 
electron with wave vector k, spin a in the lead A, V\k is the hybridization matrix element 
between conduction electron of energy E\k in a lead A with chemical potential fi\ and 
localized electron on the dot. Ed is the single particle energy of electrons in the dot. 
n-f = d^d^ is the number operator for electrons with spin up localized on the dot and U 
is the (repulsive) interaction energy between two electrons on the dot. 

The correct way to calculate the charge on the dot under non equilibrium conditions 
is to use non equilibrium Green's function (GF) of Keldysh [ [TO , RJ. In this technique the 
average charge on dot at time t is given by 

(n(t)) = £<#(*K(*)> = -i£G<(M)- (5) 

(7 cr 

where G<(t,t) is diagonal (in time indices) element of the Keldysh "lesser" GF 

G<(t,t')=i(dt(t)d a (t')). (6) 

This function is the (1, 2) matrix element of the following contour ordered Green's function 
matrix 

G(t,t') = -i(T c d{T)d + (r')) =\ Gci G< , (7) 

V G> G- c ) 

where T c is a complex time contour ordering operator. The time contour starts above 
time axis at to = — oo passes through t and t' and returns back to to — — °o below 
time axis. G c = —i(T d(t)d + (t')) is the usual time ordered (causal) GF for both times 
on the upper branch of the contour, while G e (t, t') is antitime ordered GF with both 



time arguments on the lower branch. Greater (G > ) and lesser (G < ) GFs are distribution 
functions: G > (t,t') = —i(d(t)d + (t')). Not all this functions are independent. 

It turns out that the calculation of the "lesser" GF is easy for noninteracting quan- 
tum dot i.e. for U = 0. It can be found directly in the following way. Because, the 
Keldysh contour ordered GF, possesses the same perturbation expansion as the corre- 
sponding equilibrium GF [|10j it is possible to write down the matrix Dyson equation 

G = G + G Q ZG. (8) 



The theorem due to Langreth [[Tl|] allows the transition from contour ordered to real axis 
functions and one gets || 

G< = (l + G r S r )G<(l + S a G a ) + G r S < G' a (9) 

where G r ^ denotes retarded (advanced) GF. 

For noninteracting quantum dot (U = 0) exact expression for the retarded dot GF 

reads 

G» = G5» + G5>) E» G» (10) 

with 

E>) = £ | Va^G^u;) , (11) 

\k 

where G r 0Xk (uj) — (u — £\k + iO)^ 1 is the Green's function of the lead A. From this one 
immediately finds 

E<(w) = ^V» 2 G^{ui)=i^\V Xk \ 2 8{u>-e»)f x {u>) = 

\k Xk 

= 2m^2T x (u)f x (u), (12) 

A 

where we have used symbol T x (u) to denote the average coupling of the dot to lead A; 
T x (lu) = 2vr Efc \V Xk \ 2 S(u - s xk ). Noting that G<(t, t) = J ^G < (u) and using formula (|) 
we get 

(n) = I duoG r '(u)E < (uj)G a (oj) , (13) 
which with help of Eqs. (|H]) and ( |i~2"|) we further rewrite as 

„J ETa(^) V it J 

A 



n 

a 



Note that for symmetric coupling r^(u;) = Th(lj) this formula reduces to the form (ID) with 
average distribution function f(e) as given by @. Such an expression for dot occupation 
has been used in [0 for an interacting dot. We stress, however, that it is valid for 
noninteracting and symmetrically coupled dot only. 

As one can see from the above the main difficulty in calculation of the charge on 
the dot is connected with calculation of the "lesser" self energy S < . In the noninteracting 
case the knowledge of the exact expression for self energy (|TTp allowed easy calculation of 
exact and exact expression for the charge. For interacting dot there is no way to get 
self energy exactly and one has to approximate it appropriately. One such approximation 
has been introduced by Ng |l3| and will be discussed latter on. Here we propose to use 



recently derived equation of motion method [|14] and calculate G < (uS) with the desired 
accuracy and with approximations consistent with those made in calculation of the current 
or other response functions of the system. 

To make the algebra easy we consider U — oo limit in the Hamiltonian (H) and 
use slave boson technique to handle it. We thus rewrite (^) with help of slave boson 
operators b, b + and use the commutation rules of LeGuillou and Ragoucy |L5| to evaluate 
the necessary quantum brackets. These rules treat exactly the local constraint |16|] which 
for U = oo prevents double occupancy of the dot. The procedure is simple. One rewrites 
the Hamiltonian in the form 

H SB = E(^fc - Vx)4 ka c Xka + e d Y: ftU + E Vx k (4 ka b + U + ft hcxk.) (15) 

Xka cr Xka 

and calculates the on-dot Green's function D<(u) = ((b + fa\ft^))^ using equation (c.f. 
equation (28b) of [8]) 

((A|5)) = ^<M([A5] ± )+^M(([A^]|i?))< + ^M(([A^]|5))", (16) 

with the third term of H SB taken as a interaction part Hj and lower case GFs being the 
corresponding GFs of the free Hamiltonian H (consisting of first and second terms of 
H SB ). The higher order GFs appearing at this stage have been calculated in similar way. 
In the process we have used the same kind of factorization which one uses calculating the 
on-dot retarded GF necessary to get Kondo effect. Explicitly we neglected the GFs like 
<< c\k-aCx>k'af- a b\f+b » a as they describe higher order spin correlations in the leads 



and performed the decoupling 

« ci^cu-ob+felf+b »°« f(e Xk ) « 6+/ ff |/+6 » Q 8 kk ,6 xx , . (17) 
The resulting "lesser" self-energy takes on simple form 

Z<(cu) = -2ttz £ r A (w)(l + A(o)))/ A (w) . (18) 

A 

This is our final expression for S < (c<j) to be used in calculations of the charge n. 

Now let us rederive the expression for S < (a;) using the earlier mentioned approxi- 
mation introduced by Ng [ 13j and often used in the literature |L7|. This author proposes 



to assume that Y^ioS) = AYiq(lu) and S > (u;) = ATjq where A is an unknown function 
and is noninteracting (and thus known) self-energy. For the present example it is 

given explicitly by (|i~2|). This together with the exact relations: £ > — S < = S r — S a and 



£> _ £< = - Eg leads to 



^NgM = S^-^E<(oO (19) 



or in explicit form 



£< g M = -27TZ A 2 r A(! + AH) , (20) 



Er A (^ a 

A 



which again contains the characteristic average distribution function, which for the sym- 
metric coupling reduces to f{e). Note that this formula agrees with ( |I8"D only in equilib- 
rium situation hl = /ir = \i. 

It turns out that our formula, Eq. fll8|) can be rederived by suitably generalizing the 
Ng's ansatz. To see this note that noninteracting self energy Sq (u>) above can be written 
as a sum of pieces, each of which is connected with one of the leads i.e. £q (uj) = XaXq A (u;). 
One then expects that due to locality of U term in Hamiltonian the interacting self-energy 
will also be a sum of contributions from different leads T j < (uj) = J^x^xi 00 )- Generalising 
the ansatz and writing S A (uj) = A\Y,q X independently for each A one immediately repro- 
duces result (UP . This shows that the application of equation of motion technique |L2 



is a 



correct way to to derive the "lesser" self energy. This technique, contrary to approximate 
schemes, allows mutually consistent calculations of both retarded and distribution GFs 
in Keldysh technique. 
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Figure 1: The charge on the quantum dot connected to two external electrodes in de- 
pendence of the chemical potential fiR. The parameters are the following: couplings 
Y L = T R = 17meV, chemical potential hl = and temperature T = 10~ 4 meV/kB- 

We have numerically calculated the charge n on the interacting quantum dot (with 
U = oo) coupled to leads with different chemical potentials using four formulae discussed 
in this paper. The results are presented in figures (1) and (2). Figure (1) shows the charge 
on the dot with Ed = —lmeV as a function of the position of right lead chemical potential 
fiR with hl = and for symmetric couplings = Tji(u) — T — 17meV. Note tha for 

symmetric coupling the Ng's approximation give the same results as average distribution 
function f(e) . In more general case of asymmetric coupling these two approaches lead to 
different results, but we abandoned this complication here. Solid line show results obtained 
with help of our formula (0), dashed one is calculated with help of (^), while dotted 
from (|l|) with p, given by (H). The crossing of the lines for \i L = fi R = corresponds 
to equilibrium situation. The differences between the curves get larger for increasing 
voltage V — (fiL — HR)/e applied to the system. The large (small) differences between 
various curves for negative (positive) values of is due to large (small) values of the 
on-dot density of states in the energy window between \i L and /jlr. Figure (2) shows the 
charge on the dot calculated for constant difference /i^ — [Mr = lmeV as a function of 
the parameter E4. For E& much higher or much lower than both chemical potentials the 
differences are small as expected because the density of states around /i^ and hr is small. 
The largest differences of the calculated charge are found for such values of Ed for which 
the density of states of electrons on the dot possesses the Kondo resonance. 

We conclude by stressing that the calculation of charge on the dot is a necessary step 
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Figure 2: The charge n vs Ed for constant difference fiR — = ImeV. Other parameters 
are the same as in Fig. (|TJ) . 

towards self consistent calculations of the current flowing across quantum dot devices. In 
some parameter ranges all formulae give only slightly differing results. Generally, however, 
the differences between values of n calculated by means of different formulae may be as 
large as 20-30%. The equation of motion method allows calculation of the lesser self energy 
consistent with retarded one. This is important as the current flowing in the system is, for 
proportionate couplings T L (u) = constT R (u) , expressed by the retarded Green's function 
which in turn depends on the charge n on the dot. For general couplings the consistency 
of approximations is even more important as the current depends on both retarded and 
lesser GFs and to be consistent one has to treat both self energies on equal footing. 

This work has been partially supported by the Committee for Scientific Research 
under grant 2P03B 106 17. 
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